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Abstract 

From a sample of 1172 ± 61 7r^7r^7r+ decay, we find T{D^ — > 

7r-7r+7r+)/r(Z)+ K-'K+7r+) = 0.0311±0.0018;[};^[3^^ . Using a colierent am- 
plitude analysis to fit the Dalitz plot of these decays, we find strong evidence 
that a scalar resonance of mass 478tlt ± 17 MeV/c2 and width 3241^^ ^ 21 
MeV/c^ accounts for approximately half of all decays. 
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The three-body decays of charm mesons often proceed as quasi- two-body decays 
with resonant intermediate states. In our companion paper we find that — >■ 
/o(980)7r"'" accounts for approximately half of the three-pion decay rate and — >■ 
/o(1370)7r^ accounts for more than half of what remains, clearly establishing the 
dominance of isoscalar resonances in producing the three-pion final state. In this 
paper we present a study of the singly Cabibbo-suppressed decay 7r^7r+7r+. 
We determine the ratio of decay rates T{D^ -k^-k^-k^) /T{D^ K^-k^-k^) and 
study the —>■ 7t~7i~^tt~^ Dalitz plot to determine the structure of its density 
distribution. We find that allowing an amplitude for an additional scalar state, with 
mass and width unconstrained, improves our fit substantially. The mass and width 
of the resonance found by this fit are 478^^ ± 17 MeV/c^ and 324l^g ± 21 MeV/c^. 
Referring to this tt^tt~ resonance as the (j(500), we find that (T(500)7r+ 
accounts for about half of the total decay rate. 

Other experiments have presented controversial evidence for low-mass vrvr res- 
onances in partial wave analyses P, ^ |[, with ambiguous results for the charac- 
teristics of such particles!]^, ^. Theoretically, light scalar and isoscalar resonances 
are predicted in models for spontaneous breaking of chiral symmetry such as the 
Nambu-Jona-Lasinio linear a model and its QCD extension 0]. Also, these par- 
ticles have important consequences for the quark model , for quark-gluon models 
0, for understanding low energy vrvr interactions [l^, and for understanding the 
A/ = 1/2 ruleQ. 

This study is based on the sample of 2 x 10^° events recorded in Fermilab ex- 
periment E791, in which 500 GeV/c '7r~-nucleon interactions were observed using 
an open geometry spectrometer. The final analysis makes no direct use of particle 
identification; it is solely based on tracking and vertex reconstruction capabilities. 
To reduce background, we required a 3-prong decay (secondary) vertex to be well- 
separated from the production (primary) vertex, and located well outside of the 
target foils and other solid material. The momentum vector of the D candidate had 
to point back to the primary vertex. A more detailed description of the final sample 
selection criteria and some additional details are provided in the companion paper 
m] where the resulting 7r^7r+7r+ invariant mass distribution is shown as Fig. 1. 

We have, in addition to the combinatorial background, three kinds of charm 
backgrounds: the reflection of the K^tt^tt^ decay, located below 1.85 GeV/c^ 

in this spectrum; the decay D° K^n^ plus one extra track (mostly from the 
primary vertex); and the decay chain ?7'vr^, i]' — >• p° (770)7, P°(770) vr+vr". 

The last two reflections populate the whole 7r^7r+7r+ analyzed spectrum. We use 
Monte Carlo (MC) simulations and data to determine both the shape and the size 
of each type of charm background. The combinatorial background in the vt^tt+tt^ 
invariant mass spectrum is represented by an exponential function. 
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We fit the 7r^7r+7r+ invariant mass distribution shown in Fig.l of reference |1| as 
the sum of and signals plus background. Each signal is described as the sum 
of two Gaussians with a common centroid but different widths, all these parameters 
determined by the fit. The background is represented by a function with four terms 
described above and in ^. The fit yields 1172 ± 61 D+ events and 848 ± 44 Df 
events. 

We measure the branching ratio for 7r~7r+7r"^ relative to that of —>■ 

K~TT~^7i~^. The K~7i~^7[~^ signal, selected with the same criteria used for the tt'tt^tt^ 
signal, is 34790 ± 232 events. The absolute efficiency, e, for each decay mode is ap- 
proximately 3%. From Monte Carlo studies we determine that the ratio e{D~^ 
71^71^71^) /e{D^ K^Ti^Ti^) = 1.08 ± 0.02. Note that we use the decay matrix 
element found in this analysis and an appropriate one for K^tt^-k^ to deter- 

mine the efficiencies. We also weight the MC production model to match our data. 
We find the relative branching ratio to be 

Z K-X\ = ^ °-°°18^aSS . (1) 

The first error is statistical and the second is systematic. Uncertainties in the 
7i~n~^7i~^ background shape and the levels of some contributions dominate the sys- 
tematic error. This result can be compared with the measurements reported by 
E691 [0, 0.035 ± 0.007 ± 0.003, by WA82 jg, 0.032 ± 0.011 ± 0.003, and by E687 



|T1, 0.043 ±0.003 ±0.003. 

To study the resonant structure of the decay — > 7r~7T~^TT~^, we consider the 
1686 candidates with invariant mass between 1.85 and 1.89 GeV. The integrated 
signal-to-background ratio in this range is about 2:1. Fig. |I] shows the Dalitz plot 
for these events. The horizontal and vertical axes are the squares of the tt+tt^ 
invariant masses, and the plot has been symmetrized with respect to the two 7r"^'s. 

To study the resonant structure in Fig. |^, we use MINUIT|]T6|] to extract the 
parameters. We do this by maximizing the log (Likelihood) C for several models of 
signal and background. For each model we compute C in terms of signal and back- 
ground probability distribution functions (PDF's) of the 7r^7r+7r+ invariant mass, 
M , and the Lorentz invariants Si2 = and S13 = (in our convention the odd- 
charged pion is labeled particle 1). Writing Vs and Vb for the signal and background 
PDF's, £ = n]iT[Ps + VB\r We take Vs = ^^g{M)e{s^2, ^13) | A \\ with 

N 

A = aoe*''« A + E «ne*'Mn(si2, ^13) . (2) 

n=l 

In this equation Ns is the normalization constant, e{si2i S13) is the net efficiency, 
g{M) is a Gaussian function describing the signal 7r^7r+7r+ invariant mass spectrum. 
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Figure 1: — > vr 7r"'"7r+ Dalitz plot distribution. Since we have two identical 
particles this distribution was symmetrized. 

The fit parameters are the coefficient magnitudes, a„, and the phases, 

The non-resonant amplitude, is represented by a constant. Each resonant 
amplitude, Aniji > 1), is written as a product of four terms. 

An= 'Fn^'Mn^BWn (3) 

The ffist term is form factor for the n*^ resonance We assume the form factor for 
the D decay to be one. M.n is a term which accounts for angular-momentum con- 
servation and depends on the spin J of the resonance. The final term is a relativistic 
Breit-Wigner function given by: 



m"^ — rriQ + irriQi nyn) 



with 



In Eqs. 3 and 4, m is the invariant mass of the two pions forming a spin-J resonance. 
The functions ^ F are the Blatt-Weisskopf damping factors [0: °F = 1 for spin 
particles, = l/^Jl + {rp*Y for spin 1 and = 1/ \JQ + ?>{rp*Y + {rp*Y for spin 
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2. The parameter r is the radius of the resonance (~ 3/m) and p* = p*{m) 
the momentum of decay particles at mass m, measured in the resonance rest frame, 
Po = P*("^o)5 where mo is the resonance mass. The spin part of the amphtude 
'-'Mn is defined equal to 1 for a spin-0 resonance, — 2 | ps || p2 | cos9 for spin-1 
and |(| Pa || P2 |)^(3cos^6' — 1), where pa is the 3-momentum of the unlike-charge 
pion and p2 is the 3-momentum of the other like-charge pion, both measured in 
the resonance rest frame; and 6 is the angle between pions 2 and 3. Finally, each 
resonant amplitude is Bose symmetrized: An = An[{12)3] + ^„[(13)2]. 

The background distribution is given by Vb = b{^-^)Y.i=i j^^-iisu, Sis); b{M) 
is the function describing the background distribution in the tt^tt^tt'^ spectrum, 6j 
are the relative amount of each background type and iV^. are the corresponding 
normalization constants. 

The three components of the background distribution are the combinatorial back- 
ground, assumed to be uniform before any acceptance effects, and the —>■ K~tt~^ 
and Df — > t^'tt"'" reflections. The relative background fractions are 80 ± 6%, 4 ± 1%, 
and 16 ± 6%, respectively. The shape, location, and size of the charm background 
were obtained using MC simulations and previously determined and Ds produc- 
tion rates relative to in our data sample. All parameters used for the background 
description are fixed during the fit. 

The dominance of the above three background contributions was checked in sev- 
eral tests. The analysis was repeated with more stringent selection criteria, with 
various levels of Cerenkov-counter requirements on the pions, and with varied back- 
ground levels in the fit. All test results were consistent with the quoted final results. 
In addition, Monte Carlo simulations were used to study specific charm decay chan- 
nels and a generic sample of charm decays. The latter was examined to look for 
Dalitz-plot structure in the events which passed our final selection criteria. No 
structures, other than those noted above, were significant. 

In a first model, which we will refer to as Fit 1, the signal PDF includes a 
non-resonant amplitude and amplitudes for decaying to a vr"*" and any of five 
established Tr+Tr" resonances§: p°(770), /o(980), /2(1270), /o(1370), and p°(1450). 
In the case of the /o(980) and /o(1370), we used the parameters of Ref. and 
not those of Ref. [^. The fit extracts the magnitudes and phases of each of the 
amplitudes along with the error matrix for these parameters. We calculate the decay 
fraction for each amplitude as its intensity, integrated over the Dalitz plot, divided 
by the integrated intensity of the signal's coherently summed amplitudes. The Fit 
1 results are listed in the first column of Table In this model, the non- resonant, 
the p°(1450)7r+ and the p°(770)7r+ amplitudes dominate. The qualitative features 
of this fit are similar to those reported by E691|jl3| and E687 [p!5[ . 
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Fit 1 


Fit 2 


mode 


Praction(%) 


Fraction (%) 




Magnitude 


Magnitude 




Phase 


Phase 




— 


46.3 ± 9.0 ± 2.1 
1.17 ± 0.13 ±0.06 
(205.7 d= 8.0 ± 5.2)° 


pO(770)7r+ 


20.8 ± 2.4 


33.6 ± 3.2 ± 2.2 




1 (fixed) 


1 (fixed) 




O(fixed) 


(fixed) 


NR 


38.6 ± 9.7 


7.8 ± 6.0 ± 2.7 




1 36 ± 20 


48 ± 18 ± 09 




(150.1 ± 11.5)° 


(57.3 ± 19.5 ± 5.7)° 


/o(980)7r+ 


7.4 ± 1.4 


6.2 ± 1.3 ± 0.4 




0.60 ± 0.07 


0.43 ± 0.05 ± 0.02 




(151.8 ± 16.0)° 


(165.0 ± 10.9 ± 3.4)° 


/2(1270)7r+ 


6.3 ± 1.9 


19.4 ± 2.5 ± 0.4 




0.55 ± 0.08 


0.76 ± 0.06 ± 0.03 




(102.6 ± 16.0)° 


(57.3 ± 7.5 ± 2.9)° 


/o(1370)7r+ 


10.7 ± 3.1 


2.3 ± 1.5 ± 0.8 




0.72 ± 0.12 


0.26 ± 0.09 ± 0.03 




(143.2 ± 9.7)° 


(105.4 ± 17.8 ± 0.6)° 


p°(1450)7r+ 


22.6 ± 3.7 


0.7 ± 0.7 ± 0.3 




1.04 ± 0.12 


0.14 ± 0.07 ± 0.02 



(45.8 ± 14.9)° (319.1 ± 39.0 ± 10.9)° 

Table 1: Final result with the first errors statistical and the second, in Fit 2, the 
systematic. 
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To assess the quality of the fit, we developed a fast-MC program which produces 
binned Dalitz-plot densities accounting for signal and background PDF's, including 
detector efficiency and resolution. Comparing the binned Dalitz-plot-density distri- 
bution generated by MC events using the magnitudes and phases of the amplitudes 
given in Fit 1 with that for the data, we produce the distribution for the differ- 
ence in densities and observe a concentration of a large in the low tt+tt" mass 
(mTT+TT-) region. The summed over all bins is 254 for 162 degrees of freedom (i/), 
which corresponds to a confidence level less than 10~^, assuming Gaussian errors. 
Since the two m'^+^- projections are nearly independent, we display the sum of Si2 
and Si3 in Fig. ^ for the data and for the fast-MC. 

The small value of the confidence level casts doubt on the validity of the model 
used. While the projection of the MC onto the vr+vr" mass^ axis describes the 
data in the p°(770) and /o(980) regions well, there is a discrepancy at lower mass, 
suggesting the possibility of another amplitude. 

To investigate the possibility that another tt+tt" resonance contributes an ampli- 
tude to the TT~TT^TT^ decay, we add a sixth resonant amplitude to the signal 
PDF. We allow its mass and width to float and assume a scalar angular distribu- 
tion. This fit (Fit 2) converges and finds values of 478^*^23 MeV/c^ for the mass and 
324^40 MeV/c^ for the width. We will refer to this possible state as the cr(500). The 
corresponding results, including the systematic errors , are collected in the second 
column of Table 0. 

In Fit 2, the cr(500) amplitude produces the largest decay fraction, 46%, with a 
relatively small statistical error, 9%. The non-resonant fraction, which at (39 ±10)% 
was the largest in the original fit, is now only (8 ±6)%. When we project this model 
onto the Dalitz plot, the x^/v becomes 138/162. The projection of this model 
onto the vr'^vr" invariant mass squared distribution, shown in Fig. 03, describes 
the data well, including the accumulation of events near 0.2 GeV^/c^. We tested 
the accuracy of the fit's error estimates by producing hundreds of fast-MC samples 
using the parameters of Fit 2 and then fitting the samples. The central values of 
all the parameters are reproduced accurately, and the width of the cr(500) fraction 
distribution is 0.12, slightly larger than MINUIT's estimate, 0.09. 

When comparing the two models we used the fast-MC to simulate an ensemble 
of samples for each model. For each sample we calculated Aw = — 2(ln£j — \n.C^), 
where and are the likelihood functions evaluated with the parameters from 
the fit with and without the an^ amplitude, respectively. In the data. Aw = 118; 
in the fast MC with a a, {Aw) = 108; in the fast MC with no a, (Aw) = -106. 
In both MC experiments, the rms deviation for Aw is about 20, indicating a strong 
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Figure 2: Si2 and S13 projections for data (error bars) and fast-MC (solid). The 
shaded area is the background distribution, (a) solution with the Fit 1, (b) solution 
with Fit 2. 

preference for the additional amphtude, as is indicated also by the difference in x^/z/ 
for the two models. 

We consider the systematic errors associated with the values of the fixed param- 
eters in the fit. The most important ones come from uncertainties in the background 
model (the background shape, composition, and level), in particular the D'^ — > ryV"'" 
reflection which populates the same region as the D'^ — > p^Tr"*" component. We also 
account for the uncertainties in the parameters describing the acceptance function. 

To better understand our data, we also fit it with vector, tensor, and toy models 
for the sixth (sigma) amplitude, allowing the masses, widths, and relative amplitudes 
to float freely. The vector and tensor models test the angular distribution of the sig- 
nal. The toy model tests the phase variation expected of a Breit-Wigner amplitude 
by substituting a constant relative phase. The vector resonance model converges to 
poorly defined values of the mass and width: 805 ± 194 MeV/c^ and 1438 ± 903 
McV/c^; the tensor model to more poorly defined values: 2350 ± 683 McV/c^ and 
-690 ± 1033; and the toy model to 434 ± 11 MeV/c^ and 267 ± 37 MeV/c^. As 
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a test of the models, we again project the vector, tensor, and toy models onto the 
Dalitz plot and obtain ^jv = 188/162, 148/162, and 152/162, respectively. For 
these models, we also find Aw = 66, 13, and 15 where MC experiments predict 
(Aw) ~ 64, 56, and 38 when the data is generated with the scalar parameters and 
the negatives of those values when the MC data is generated according to the vec- 
tor, tensor, and toy model parameters. The rms widths of the MC distributions are 
15, 15, and 11 units respectively. These statistical tests strongly exclude the vector 
model. They clearly prefer the scalar model to the tensor and toy models. Note 
that the central value for the tensor mass is well above threshold for decay and 
the negative width is an indication that no physically meaningful tensor resonance 
fits the data. In the toy model, the extra amplitude interferes strongly with a large 
non-resonant amplitude, leading to an unphysically large sum of resonant fractions. 

In summary, from 1172 ± 61 '7r~7r"'"7r"'" we have measured r{D~^ — >■ 

7r-7r+7r+)/r(D+ ^ K'n+TT+) = 0.0311 ± 0.0018l|];g[]^^. In an amplitude analy- 
sis of a sample with S:B ^ 2:1 we find strong evidence that a scalar resonance with 
mass 478^23 ± 17 MeV/c^ and width 324l4o ± 21 MeV/c^ produces a decay fraction 
~ 50%. Alternative explanations of the data fail to describe it as well. The promi- 
nence of an amplitude for an isoscalar plus a vr"*" in this decay accords well with our 
observation that the amplitude for an isoscalar plus a pion produce a large majority 
of all TT^TT+TT"'" decays[|ll. 
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